Arbitrary rotation and entanglement of flux SQUID qubits 
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We propose a new approach for the arbitrary rotation of a three-level SQUID qubit and describe 
a new strategy for the creation of coherence transfer and entangled states between two three-level 
SQUID qubits. The former is succeeded by exploring the coupled-uncoupled states of the system 
when irradiated with two microwave pulses, and the latter is succeeded by placing the SQUID qubits 
into a microwave cavity and used adiabatic passage methods for their manipulation. 
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I. INTRODUCTION 

It has been realized over the last few years that solid 
state systems that make use of the Josephson effect 
could play an important role in the area of quantum 
computation^.. A series of successfully performed in- 
teresting experimentsSiSi^iSiSiLSiSiiSiiLiS have confirmed 
the applicability of these systems. Among superconduct- 
ing systems performing quantum computations, empha- 
sis has been given to the study of schemes based on mag- 
netic flux states in superconducting quantum interference 
devices (SQUID)iii3ii^.i5,i6,i7,is,i9,20^ 

Zhou et alm^ have recently proposed a three-level A- 
type rf-SQUID qubit. Here, the states of the qubit are 
the two lower flux states |0) and |1) of the SQUID sys- 
tem that reside in two distinct potential valleys, see flgure 
^ The manipulation of the qubit is done with two mi- 
crowave fields that couple the lower states to an upper 
state |e) in a A configuration. As the coupling matrix 
elements corresponding to the transitions |0) ^ |e) and 
|1) ^ |e) are larger than that of the |0) ^ |1) transition, 
the three-level SQUID qubit has been shown to be more 
favorable than the conventional two-level SQUID qubit. 
Amin et al.^^ have latter shown that the approach of 
Zhou et al. was incomplete and have proposed a more 
general method for producing an arbitrary qubit rotation 
using the three-level SQUID qubit. In addition, more re- 
cently Yang and Han^^ have shown that large detuning of 
the driving fields from the upper state could be favorable 
for implementing single qubit rotation in the three-level 
SQUID qubit. Finally, Yang et alM^Si have proposed 
two different strategies using three-level SQUID qubits 
in a microwave cavity for entanglement, logical quantum 
gates and quantum information transfer. 

In this article our goal is two- fold. First, we propose a 
new approach for the arbitrary rotation of a three-level 
SQUID qubit. The basic idea of the rotation procedure 
described here is related to that studied by one of the 
authorsSZ for atoms driven by optical fields. However, 
there the rotation is performed by means of an adia- 
batic process, that requires the involvement of an ad- 
ditional long-lived auxiliary level and an extra coupling 
field. Here, we use only two coupling fields and present a 



rotation method based on Rabi oscillations. As we are in 
the microwave domain, the sufficiently precise control of 
the pulse area is possible, unlike in the optical domain. 
Second, we apply a new strategy for the realization of co- 
herence transfer and for the creation of entangled states 
between two three-level SQUID qubits. This is succeeded 
by placing the qubits into a microwave cavity and by us- 
ing adiabatic passage methods for their manipulation. 
Both of our approaches are fundamentally different from 
the ones that have been proposed so far in the litera- 
ture as they are based on the exploration of the coupled 
and uncoupled states of the system. Moreover, the latter 
method is an adiabatic method, and such methods are 
robust with respect to fluctuations of several experimen- 
tal parameters'^. We will discuss this issue for our case 
below. 

The article is organized as follows. In the next sec- 
tion we briefly summarize the model for the three-level 
A-type SQUID qubit-Ji and present our method for arbi- 
trary qubit rotation. The general method is also special- 
ized in two simple cases. Then, in section HI we study 
the interaction of two three-level SQUID qubits in a mi- 
crowave cavity and show that using adiabatic methods 
both robust quantum information transfer and entangle- 
ment are possible. Finally, we summarize our results in 
section IV. 



II. ROTATION OF A THREE-LEVEL SQUID 
QUBIT 

The model qubit consists of an rf-SQUID which inter- 
acts with two microwave fields. The rf-SQUID is made of 
a superconducting ring interrupted by a Josephson tunnel 
junction. The system's generalized coordinate is the total 
magnetic fiux in the ring $, subjected to the potential24 
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Here, L is the ring inductance, is an externally ap- 
plied magnetic flux to the SQUID, Ej = Ic^o/2tt is the 
maximum value of the Josephson energy, with Ic being 
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the critical current of the junction, and $o = h/2e is the 
flux quantum. The Hamiltonian of the system can be 
written as 



(2) 



where Q = —ihd/d^ is the charge on the junction capac- 
itance C. The flux $ and the charge Q are canonicaUy 
conjugate operators satisfying the commutation relation 
[$, Q] = ih. A typical potential of this form is shown in 
Fig.d 

We will flrst discuss the case of single SQUID qubit 
rotation. To achieve this the SQUID is driven by two mi- 
crowave pulses. The microwave pulses are considered as 
linearly polarized electromagnetic fields with their mag- 
netic fields perpendicular to the plane of the SQUID ring. 
We take the angular frequencies of the two microwave 
fields (Do and wi to be near resonant with the transitions 
|0) ^ |e) and |1) ^ |e) respectively, where |0) and |1) 
are the states of the SQUID qubit and |e) is an excited 
state, as is shown in figure ^ The resonant approxima- 
tion can then be used to describe the dynamics of the 
system, i.e. we assume that only the states |0), and 
|e) contribute to the dynamics of the system. This ap- 
proximation has already been used successfully in several 
recent articles and the system has been termed three-level 
A-type SQUID qubitiii2ii2*22i2i. In order to analyze our 
system further we use the interaction picture and apply 
the rotating wave approximation^^i^. Then, the Hamil- 
tonian of the system can be expressed as 



H 



(3) 



where A,- 



LUj, with j = 0, 1 is the microwave 



field detuning from resonance with the |i)<->|e) transition, 
where hujj denotes the energy of the jth stationary state 
of the SQUID. 

We assume that Ao = Ai such that the system is 
at two-photon resonance. Then, in the rotating wave 
picture^, obtained by applying a unitary transformation 
to the Hamiltonian of Eq. the Hamiltonian of the 
system is given by 



H' ^hA\e){e\ + -[no\0){e\ 



17i|l)(e| +H.C.; 



(4) 



We require that the light pulses share the same time- 
dependence but their peak amplitudes can be differ- 
ent, and there can be a phase-difference between them. 
Hence, the Rabi frequencies in the Hamiltonian of Eq. 
read 

^oit) ^n{t)cos(p, Qiit) ^fl{t)e"^smip, (5) 

where rj and f are fixed angles. These two pulses define 
a coupled state |C) 



and an uncoupled state 

\NC) = -sin(^|0) +e*''cos(^|l) 



(7) 



with respect to the microwave pulsesSi. In this basis the 
initial state \4>i) of the SQUID under consideration is 
given by 

= {NC\^P,)\NC) + {C\^P,)\C) , (8) 
and the Hamiltonian Eq. ^ reads 

H' ^hA\e){e\+'^[n\C){e\+ll.c.] . (9) 

This Hamiltonian is that of a two-level system. There are 
several analytically solvable models for two-level systems 
with pulsed excitation. The transfer matrix of the general 
solution can be parameterized as 
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in the basis {\NC), |C), |e)}. The columns of the matrix 
tj{tf,ti) correspond to the state vector of the system at 
time tf if it was initially in states \NC),\C), and |e), 
respectively. Here, we need such a model that after the 
pulse has passed the coupled state |C) acquires a phase 
shift ~S and the excited state |e) is not populated. Hence 
the parameters of the transfer matrix U{tf,ti) should be 
given by 



/3 = 0. 



(11) 



Therefore, at the end of the process the state of the 
SQUID is given by 

= {NC\i^,)\NC)+e-''{C\i,,)\C). (12) 

By inserting the explicit form of the scalar products 
{NC\ijji) and {C\ijji), and the definitions Eqs. ®, 
into Eq. 112|) we obtain 



(13) 



\C) =cosv3|0)-He"'sin(/j|l) 



(6) 



where n = (sin 2ip cos ry, sin 2ip sin 77, cos 2lp) . Apart from a 
global phase —5/2 the states and 1-0/) are connected 
through the rotation Rn{S). The rotation Rn{S) is an el- 
ement of the SU(2) group and describes a rotation about 
the axis n, through the angle S. If the qubit is isolated, 
then the global phase —5/2 is unimportant. If the qubit 
is part of a larger system, e.g. there are several qubits 
which form a quantum computer, then the global phase is 
clearly relevant, however, it may be incorporated into the 
algorithm being implemented on the quantum computer. 

The simplest model for a two-level system that can be 
used to realize the dynamics described above is the Rabi 
model with rectangular pulse shape and constant detun- 
ing. Another possibility is the Rosen-Zener model22i2S 
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with hyperbolic-secant pulse shape and with constant de- 
tuning as well. In the case of the Rosen-Zener model the 
phases are given by rather involved formulae and we will 
not discuss it here. A further, and rather general model, 
is the one obtained under pulsed excitation in the case 
of far off-resonant Raman coupling, i.e. the case that 
A ^ |J7(t)|/2. Below, we describe briefly the Rabi model 
and the off-resonant Raman model. In the Rabi model 
the elements of the transfer matrix (|10|l are given by 



cos I if2T 



A . /1~ 
i— sin —flT 

n V2 



-i— sin -ilT I e 

n V2 



iAT/2 



e-'^^/2(14a) 
(14b) 



where Q = vfP + A^. For pulse area flT = 2-Km, with 
m being an integer, the transition amplitude (3 is zero, 
as is required in Eq. Hll|l . Hence the rotation angle 5 is 
given by 



(15) 



We note that the actual physical system is described by 
the Hamiltonian 10} . The two Rabi frequencies Slo and 
ill are derived from the same pulse according to Eq. ||SJ), 
which is described by the Rabi frequency Q, detuning A, 
and duration T calculated above. 

In the off-resonant Raman modeliS state |e) is elimi- 
nated adiabaticall\i2Si^ and 



a « exp I — m)\^dt 
/3 « 0. 

Therefore, the rotation angle reads 



1 

'4A 



\n{t)\^dt. 



(16) 
(17) 

(18) 



We have performed simulations for a realistic SQUID 
system. We used the same parameters for the SQUID as 
in the work of Zhou et al^, i.e. L = 100 pH, C = 40 fF, 
Ic = 3.95fj.A and = -0.501$o- For this SQUID the 
dissipation time could exceed 1 ^seoi. The qubit rotation 
time was found to be dependent on the model of interac- 
tion that we used and varied from sub-nanosecond times 
to about 30 nsec for moderate Rabi frequencies values of 
maximum strength of 1-5 GHz. An example of inversion 
from state |0) to state |1) using the Rabi model is shown 
in Fig. [3 We note that shorter qubit rotation times can 
be achieved by increasing the microwave field intensities. 
Of course, this cannot happen arbitrarily as after a cer- 
tain limit the rotating wave approximation will not be 
appropriate for describing the system and the effect of 
other states, that have been omitted here, will have to 
be taken into account. 



III. QUANTUM INFORMATION TRANSFER 
AND CREATION OF ENTANGLED STATES 

We will now present a strategy for achieving entangle- 
ment and also information (coherence) transfer between 
two A-type SQUID qubits. We place both SQUIDs (we 
denote them by A, B) in a microwave cavity, see Fig. 



13 and assume that the transitions |l)j 



\e')j, with 



J 



A, B and where 



are excited states different from 



those used for the single qubit rotation, are coupled to 
the same cavity mode, with coupling constant g, which 
is assumed real. The other transitions |0)j |e')j-, with 
j = A,B are coupled with external laser fields, with Rabi 
frequency Qj, with j — A,B, such that each microwave 
field addresses individually only one SQUID. We assume 
again that the cavity field and the external microwave 
fields are at two photon resonance. The Hamiltonian for 
the two SQUID system in the rotating wave picture and 
in the rotating wave approximation is given hyiS, 



(19) 



Hab = A'\e'),,{e'\ 

j=A,B 

+ f E {^j\0)n{e'\+g\l)jj{e'\b^+n.c 

3=A,B 



where W is the creation operator of a cavity photon. 
We are interested in the uncoupled states of the two 
SQUID system, which involves the vacuum cavity state 
|0)c- These are 



(I7^i(i)g|l,0) + f7|j(t).g|0,l))|0)e 



-n^mBim^A) 

|1,1,0), 



(20a) 
(20b) 



where N = 1/ v/pIFTp^P)?TpIP^W- These 
states are uncoupled because 



HAB\i^'' 



0, 



1 



I. II. 



(21) 



The state of Eq. (j20b|l is constant, it doesn't change with 
time. The other uncoupled state Eq. H20a(l may vary with 
time as the Rabi frequencies 0,A{t) and ^B{t) can be time 
dependent. 
The subset 

Wo = {|0,l,0),|e',l,0),|l,l,l),|l,e',0),|l,0,0)} (22) 

of the basis states of the system is closed in the sense that 
the matrix elements of the Hamiltonian Eq. (|19|l between 
any other basis state and a state taken from this subset is 
zero. Here, \a,b,n) = | a) ^| 6) B|n)c where |n)c denotes the 
state of the photons in the cavity. Since the uncoupled 
state Eq. (|20a|l is a superposition of the basis states from 
the subset Eq. (|22|) , then the dynamics of this uncoupled 
state can be obtain in the subset TYq- In the basis states 
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of Ho the Hamiltonian (|19|l is given by 



H'ab - 



r^A 

n\ 2A' 5 

5 5 

5 2A' n*g 

flB 



(23) 



Let us consider first information transfer between 
SQUIDs A and B^. If the system is prepared initially 
in a state of the form 



li^iU)) = (co|OU + ci|1)a)|1)s|0), 



(24) 



with Co and ci being arbitrary coefficients satisfying 



1, then it is said that the state of the qubit 



A is transferred in qubit B if after some interaction the 
state of the system is given by 



|V'(i/)) = |1)a(co|0)b+ci|1)s)|0). 



(25) 



How can we realize this process in the two SQUID sys- 
tem? We observe that the uncoupled state Eq. H2()a|l 
is equal to the state |0, 1,0), for fig V,a- Hence let 
us chose the external pulses A and B such that at early 
times around ti the Rabi frequencies satisfy the condition 
^3(^1) ^ ^A{ti)- Then the initial state of the system 
Eq. 124|l can be expressed as 



|V'(^^)) =Co|V/) +C1IV' 



(26) 



In the opposite limit when ^A, the micoupled 

state coincides with the state |1, 0, 0). Let us choose 
the time dependence of the pulses A and B such, that 
at early times we have flsiti) ^ ^aIU), then after 
a smooth variation, at late times we have ^Isitf) ^ 
^A{tf). We also require that the two pulses overlap. 
These conditions ensure that the uncoupled state is 
well-defined throughout the whole time evolution and the 
state of the system follows adiabatically this uncoupled 
staters. The pulse sequence described here resembles 
that of the stimulated Raman adiabatic passage (STI- 
RAP), for reviews see^. Hence, in the above adiabatic 
process the state co\ijj^) -|-ci|V'^^) will evolve smoothly to 
the final state given by Eq. (|25|l . and the required infor- 
mation transfer between SQUIDs A and B will be real- 
ized in this way. The adiabaticity conditions are given 
by 



|(^fe|^')|«k/c| 



k = i. 



(27) 



where \^k) is an instantaneous eigenstate of the Hamil- 
tonian H23|l . and Ek is the corresponding non-zero eigen- 
value 



(28) 



The fifth eigenstate is \ip^) belonging to the eigenvalue 
zero. The eigenvalues £k are given by 



£k 



2 ± ^ V4A2 + ± 2^{\nA\' - |r!sP)2 + 4.g4 , 

(29) 



with rt^ = 2g^ + I^IaI"^ + \^b\'^ , and the eigenvectors read 



1 



g^^A 

-g{\nA\^+4ek{A-ek)) 
-2£fe(.92 + |f7^|2 + 4e,(A-e,)) 

-nB{g^ + \nA\^ + iekiA-ek)) 



(30) 



By inserting Eqs. ^U^, and ^ into Eq. ^ the 

fulfillment of adiabaticity can be verified for the chosen 
external Rabi frequencies, cavity coupling strength, and 
detuning: similarly to the standard three-level STIRAP, 
large pulse areas (fi^i.sT, gT ^ 1, with T being the char- 
acteristic length of the pulses), and smooth, slowly vary- 
ing pulse envelopes are required. 

We now address the question of creating entanglement 
between the two SQUID systems. The entangled state is 
defined by 



|V(^,0> =(cos0|l,O) 



e-'^ sin 6*10, 1))|0), 



(31) 



where 9 and ^ are fixed angles. In order to ob- 
tain this state, we apply the method of fractional 
STIRAP34.35.36,37^ ^j^ogg ^^le initial state of the sys- 
tem to be IV'(^i)) — lOjljO) and the order of the pulses 
as before, i.e. at early times {t ^ ti) flBit) ^ ^A{t), so 
that the uncoupled state coincides with the initial 
state of the system. If the evolution is adiabatic and the 
two pulses are switched off such that flA{t) = ^{t) cos 9 
and ris(i) = il{t) sm6e^^ as i ^ tf, then the entangled 
state Eq. H31|) is created. Moreover, by means of the 
inversion of the state of SQUID A, the entangled state 



\^{d,0) = (cos0|O,O) 



e-'« sin 6*11, l))|0)c 



(32) 



can also be generated. 

Our information transfer and entanglement scheme 
possesses the merits of both SQUID cavity QED quan- 
tum computing schemes2P and adiabatic evolution trans- 
fer schemes^. The major advantages of SQUID cavity 
QED schemes are the simplicity of the coupling between 
the two SQUIDs via the cavity mode, the protection of 
SQUIDs from the interaction with the environment and 
thus the reduction of decoherence, and the technical sim- 
plicity of placing SQUID in cavities. The main advantage 
of adiabatic evolution transfer methods, such as STIRAP, 
is the robustness of the method for moderate fiuctuations 
of the microwave pulse parameters. In addition, as long 
as the evolution remains adiabatic the excited states \e')j , 
with j = A, B are minimally populated, while the cavity 
mode is only populated in the transient regime. The lat- 
ter can be greatly suppressed, or completely avoided, by 
keeping the coupling g greater than the Rabi frequencies 
rij, with j — A,B. Therefore dissipation and decoher- 
ence can be reduced using this scheme. 

We have performed simulations for the SQUID pre- 
sented in section II. For coupling strengths for the Rabi 
frequencies and the cavity coupling coefficient of the or- 
der of 1-5 GHz information transfer takes from 10-30 
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nsec while the creation of entangled states needs slightly 
larger times up to 50 nsec. A typical example for switch- 
ing from state |0, 1,0) to state |1,0, 0) is shown in Fig. 
01 We also present the creation of the entangled state 
1/V2(|0,1,0) + |1,0,0)) in Fig.E 

IV. CONCLUSION 

In summary, we have presented schemes for basic state- 
manipulations of a single rf-SQUID and two rf-SQUID 
qubits. In the first part of the paper we have consid- 
ered a scheme for arbitrary rotation of a three-level rf- 
SQUID qubit. The rotation is performed by irradiating 
the SQUID with two microwave pulses that share the 
same time dependence. These pulses define a coupled and 
an uncoupled state out of the two SQUID states forming 
the qubit, and the rotation of the qubit results from the 
phase-shift on the coupled state caused by the applied 
microwave pulses. The main advantage of our proposed 
scheme is that it can be implemented for a large variety of 
pulse shapes and detunings, according to the capabilities 
of the experimentalists. 

In the second part of the paper we have proposed 
schemes for information transfer and creation of en- 
tangled states between two rf-SQUIDs. These pro- 



cesses are mediated by a microwave cavity, so that the 
SQUIDs communicate via photon exchange through a 
cavity mode. We have applied an adiabatic population 
transfer scheme, therefore our method is robust with re- 
spect to the moderate fluctuations of the experimental 
parameters. For proper choice of the external Rabi fre- 
quencies the excitation of the cavity field can be almost 
completely suppressed, hence the decoherence due to the 
imperfection of the cavity can be minimized. 

In closing, we note that our schemes are quite general 
and can be applied to other areas of quantum computa- 
tion where three-level qubits are used such as, for exam- 
ple, in atoms and trapped ions in cavities^S*^. 
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FIG. 2: The time evolution of the populations in states |0) 
(solid curve), |1) (dashed curve) and |e) (dot-dashed curve) for 
parameters f2 = 2 GHz, ip = 57r/4, 77 = tt, m = 2, A = —2/y/S 
GHz and 5 = tt, using the Rabi model. 



FIG. 1: The potential energy, the stationary states and the 

coupling configuration of the rf-SQUID system. The states 
|0), |1), and |e) are coupled my means of two microwave pulses 
with Rabi frequencies f2o,i in a A configuration. 




FIG. 3: A schematic representation of two rf-SQUIDs in a 
microwave cavity. The two SQUIDs are addressed individu- 
ally by microwave pulses. They are also coupled to the cavity 
field with coupling strength g. 
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FIG. 4: (a) The Rabi frequencies (dashed curve) and 

(solid curve) for the case that ftA{t) = Qe"'*"^'*^^''^^ , 
flB{f} = nf:'"'*"^s''/^p for parameters A = -2 GHz, g ^ 3 
GHz, TA = 23 ns, tb = 17 ns, Tp = 6.5 ns. (b) The time 
evolution of the populations in states |0, 1,0) (solid curve), 
|1, 0, 0) (dashed curve) and |1, 1, 1) (dot-dashed curve) for the 
above parameters. 
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FIG. 5: (a) The Rabi frequencies |nA(t)| 

(dashed curve) and (solid curve) for the 

case that iiAit) = Hsinee"**"^-*'''''^?, nsit) = 
n ^-(^-'Bflrl j^^Q^Q^-{t-TA?l-rl'^ foj. parameters = -2 

GHz, g = 3 GHz, ta = 38.5 ns, tb = 25 ns, Tp = 10 ns, 
9 = 7r/4. (b) The time evolution of the populations in states 
|0, 1,0) (solid curve), 1 1,0,0) (dashed curve) and 1 1,1,1) 
(dot-dashed curve) for the above parameters. 



